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ABSTRACT 

(-H ' We study the formation and dynamical evolution of clusters with multiple stellar 

generations. Observational studies have found that some globular clusters host a pop- 
I ' ulation of second generation (SG) stars which show chemical anomalies and must have 

Q I formed from gas containing matter processed in the envelopes of first generation (FG) 

cluster stars. We study the SG formation process by means of ID hydrodynamical 
simulations, starting from a FG already in place and assuming that the SG is formed 
by the gas ejected by the Asymptotic Giant Branch (AGB) stars. This gas collects 
in a cooling flow into the cluster core, where it forms SG stars. The SG subsystem 
emerging from this process is initially strongly concentrated in the cluster innermost 
regions and its structural properties are largely independent of the FG initial proper- 
ties. We also present the results of a model in which pristine gas contributes to the SG 
formation. In this model a very helium-rich SG population and one with a moderate 
helium enrichment form; the resulting SG bimodal helium distribution resembles that 
observed for SG stars in NGC 2808. 

By means of iV-body simulations, we then study the two-population cluster dy- 
namical evolution and mass loss. In our simulations, a large fraction of FG stars are 
lost early in the cluster evolution due to the expansion and stripping of the cluster 
outer layers resulting from early mass loss associated with FG SN ejecta. The SG pop- 
ulation, initially concentrated in the innermost cluster regions, is largely unscathed by 
this early mass loss, and this early evolution leads to values of the number ratio of SG 
to FG stars consistent with observations. We also demonstrate possible evolutionary 
■ routes leading to the loss of most of the FG population, leaving an SG-dominated 

cluster. As the cluster evolves and the two populations mix, the local ratio of SG to 
FG stars, initially a decreasing function of radius, tends to a constant value in the 
inner parts of the cluster. Until mixing is complete, the radial profile of this number 
ratio is characterized by a flat inner part and a declining portion in the outer cluster 
regions. 

Key words: globular clusters: general, hydrodynamics, N-body simulations, stars: 
chemically peculiar 
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1 INTRODUCTION 

The common paradigm that Globular Clusters (GCs) are ex- 
amples of "Simple Stellar Populations" (SSP), assemblies of 

* E-mail:annibaIe. dercole@oabo.inaf.it, vesperin@physics.drexel.edu, 
dantona@mporzio.astro.it, steve@physics.drexel.edu, 
simone . recchi@univie . ac . at 



stars born all at the same time and sharing identical chem- 
ical composition, had until a few years ago only one notice- 
able exception: uj Cen, the most massive cluster, whose stars 
show a large spread in metallicity. Today, the situation is 
much more complex, and it is now clear that practically none 
of the GCs so far examined in detail is an SSP. Even when all 
the stars in a single GC share the same Fe abundance, star- 
to-star differences in the abundances of lighter elements are 
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found. These differences concern light-element abundances: 
in a relevant fraction of cluster stars, they look like the re- 
sult of nuclear processing through proton capture reactions 
at high temperature (CN, ON, NeNa and MgAl cycles). In 
particular, O-Na and Mg-Al anticorrelations are found for 
stars at different stages of their evoluti on, from the red gian t 
branch to the main sequence (MS) (e.g. lGratton et al.ll2004) . 
Halo field stars with the sam e metallicity do not show such 
huge variations (Kraft"l994). This leads us to attribute the 
anomalous composition to the fact that these stars formed 
in a dense GC environment, and that nuclear processing in 
a first stellar generation (hereafter FG), has been followed 
by a subsequent episode of star formation, creating a second 
stellar generation (hereafter SG) from matter including this 
processed gas. 

An additional ingredient to the complexity of the new 
scenario has been added by recognizing that there are at 
least two clusters, uj Gen and NGC 2808, in which the 
MS morphology shows that a non negligible fraction of the 
stars (-25% in uj Gen , ~15% in NGC 2808) populate a 
"blue M S" which can only be interpreted as a very helium- 
rich MS llBed in ct al."20o3: lNorrisll2004l : iPiotto et al]|2005l : 
iD'Antona et a l. 20051j). The extreme helium population 
must also be very homogeneous, both in helium and metal 
content, since in both clusters the blue MS is w e ll det ached 
from the rest of the main sequence. iPiotto et all (|2007l ) have 
shown that the MS of NGG 2808 is actually made up of 
three different sequences: a normal MS, a population with 
somewhat enriched helium (Y~ 0.30) and a very helium-rich 
one (Y--0.35-0.40). In uj Cen the "normal" MS is actually 
made up of sequences with different metal content, as also 
shown by the compl exity of the subgiant branch(es) (e.g. 
IVillanova et al.ll2007h . In addition to u Gen and NGG 2808, 
the two massive clusters NGG 6441 and NGG 6388 appar- 
ently also harbor a very high helium population, including 
~ 10% and ~ 20% of stars respectively. In this case the 
evidence comes from the high- Teg extension of the Hori- 
zonta l Branch (HB) JCaloi fc D'Antonall2007l : IBusso et all 
I2OO7I : lYoon et ahl I2OO8I ), a very p eculiar occurrenc e at the 
high metallicity of these clusters (|Rich et al.lll997i ). and it 
is not clear whether the very high helium population is as 
homogeneous as it is in the previous cases. 

Observationally, it appears that these extreme Y se- 
quences are prese nt only in the most massive clusters 
llPiotto et al. !l2007'). A moderately helium rich population 
(0.26 < Y < 0.30), on the other hand, is present in most 
GGs, largely independent of the total cluster mass. The ob- 
servational spectroscopic evidence resides mainly in the ex- 
istence of similar Na-O anticorrelation — also among the 
unevolved stars — in al l the clusters so far examined (see, 
e.g., the summary by iGarretta et al.l I2OO6I ) although the 
most extreme Na-O anomalies are present only in clusters 
that also show extreme high-Tcff extensions of the HB (see 
IGarretta et al . 2007). These results, based generally on spec- 
troscopy of a few stars, have been recently confirmed by the 
analysis of larger samples of 30-100 stars per cluster (Gar- 
retta et al., 2008) The indications coming from the analysis 
of the GNO anom alies are similar (for a summary, see e.g. 
iGohen et al]|2005l ). 

Photometrically, this non-extreme population can not 
be easily isolated in the MS photometry, mainly because 
the moderately helium rich isochrones merge with the stan- 



dard isochrones in the turnoff region l|D'Antona et al.ll2002l : 
ISalaris et aLll2006l ). but also because, in ge neral, these SG 
stars do not have an uni que value of Y (e.g. iLee et al ] |2005l : 
ICaloi &: D'Antoiiall2007l ). Nevertheless, these stars show up 
quite prominently in the HB morphology, which ampli- 
fies the mass differences (due to t he initial helium con- 
tent) among stars of th e same age (|D'Antona et al.ll2002l : 
iD'Antona fc Caloilliooih Fl It is importan t to notice that, ac- 
cording to the detailed H B modelling by lGaloi fc D'Antonal 
(2007.) (see also ID'Antona fc Calol 200a ). also in NGG 6441 
and NGC 6388 the total SG population includes ~ 60% of 
the total number of stars, and is therefore much larger than 
just the very high helium fraction population. 

The interpretation of the HB morphologies in terms of 
helium enrichment, coupled with the spectroscopic informa- 
tion indicating that the helium-rich samples probably coin- 
cide with t he sample of stars showing the Na-O anomalies 
(e.g. ICarre tta et al. 2006), lead us to conclude that, in the 
clusters so far examin ed, the percentage of sta rs of the SG 
is generally -40-60% (|D'Antona fc Galoill200d ). 

For the clusters having mainly blue HBs, there are even 
indicatio ns that th ey may contain only — or mainly — SG 
stars (D'Antona fc C aloi 2008). In an analysis of the relative 
mag nitude locations of the t urnoff, red giant "bump" and 
HB. lGaloi fc D'Antonal l|2005l ) have provocatively suggested 
this possibility for M 13, the cluster showing th e broadest 
range of chemical anomalies (|Sneden et al.|[2004h . In addi- 
tion, some clusters in which the tightness of the evolution- 
ary sequences in t he HR diagram (e . g. the relatively sm all 
cluster NGG 6397, iKing et al.|[l998l : iRicher et al.ll2006l ) or 
the short extent of the HB (e.g. the massive GC M 53, 
ID'Antona fc Caioil |2008') indicate that they indeed repre- 
sent SSPs, might be dominated by a homogeneous SG. For 
example, in NGC 6397 o nly three out of fourt een scarcely 
evolved stars examined bv lCarretta et aL I (|2005l ) are normal 
in their nitrogen content. In the others, [N/Fe]~ 1 — 1.7, 
showing GN and ON processing in stars that should not 
show it. The b lue HB of M 53 can eas ily be described as 
uniparametric (|D'Antona fc Caioil l2008h . but its integrated 
spectrum shows a strong nitrogen enhancement with respect 
to field halo stars (|Li fc B urstcin 2003). In conclusion, some 
"normal" blue-HB GGs (with typical current masses of a 
few xIO^Mq ) should be considered self-enriched, and any 
model for the formation and evolution of multiple popu- 
lation clusters must be able to explain the observed large 
SG population and the cluster-to-cluster variations in the 
SG/FG number ratio. 

Both massive Asymptotic Giant Branch (A GB) stars 
ICottrell fc Da Costal Il98ll : IVentura et ai] I2OOII) or "fast 
rotating" massive star s fFRMS) (|Maeder fc Mevnetll2006l : 
iDecressin et al.l l2007al ) have been proposed as possible 
progenitors of the SG. The most critical point, common 
to both scenarios, is the following: for a standard Ini- 
tial Mass Function (IMF) of the FG, such as a Salpeter 
or Kroupa IMF, the gas ejected by the massive stars 
during the phase of H-burning, or the gas contained in 
the envelopes of massive AGBs (ID'Antona fc Galoil liooi : 



^ The age difference between FG and SG is in any case so small 
with respect to the total GC age that it is irrelevant for the turnoff 
evolution. 
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iBekki fc Norrisll2006l ), is too scarce to form a large SG pop- 
ulation. Two solutions to this conundrum have been sug- 
gested: either the system had a standard IMF but it was 
initially much more massive — at least a factor 10 more mas- 



than today's GC ()D'Antona & Caloi 


I2OO4I: iBekki & Norrid 


20061: Prantzos & Charbonnell l2006l: 


D'Antona et al. l2007l: 


Decressin et al. 2007bl), or the system evolves at constant 



l2004l : lDecressin et al.ll2007bi . see also lDowning fc Sillsl (120071 ) 



for a dynamical study) and contains a significant number of 
FG polluters. As we will discuss in Sect. 14.31 a large initial 
mass is necessary also in this second scenario if a cluster 
must have today an approximately equal number of long- 
lived SG and FG stars (where we define long-lived stars as 
those having initial masses in the range 0.1 Mq-0.8Mq). 

We reconsider here the problem of GC formation to ex- 
plain not only the most peculiar very helium-rich population 
(15-20% in three GCs, 25% in uj Gen ), but the whole SG 
population (about 40-60% of the current GC stellar content, 
and possibly even more in the "blue HB" clusters). 

In this work, we follow the idea that the intermediate- 
mass stars of the FG, starting at SM© and extending at most 
down to ~4-5M0, lose by stellar winds and the ejection of 
planetary nebulae matter processed by the hot CNO cycle 
via Hot Bottom Burning (Na- and Al-rich, O- and Mg-poor) 
during their AGB evolution, and that an SG forms either 
from the pure AGB ejecta or from ejecta mixed with pristine 
gas. While previous work examined the reliability of this 
model from a chemical point of view, here we focus mainly 
on its hydrodynamical and dynamical consistency. 

In Sect. |2] we first briefly summarize the state of the art 
of chemical modelling of the SG, we discuss the reasons why 
we focus our attention on the scenario in which the AGB 
stars are the polluters and emphasize the possible role of 
super- AGB stars. 

In Sect. |3] and [4l we address this problem by means of 
ID hydro simulations. We show that in the massive AGB 
epoch, the slow winds collect in a cooling flow driving AGB 
ejecta to the center of the cluster, where they form SG stars. 
Simulations exploring the dependence of our results on the 
cluster mass, IMF, initial concentration as well as the role 
of SNe la and other extra energy sources are also presented 
in this section. 

In Sect.[S]we explore a model in which the pristine mat- 
ter is partially re-accreted after the SN II explosion stage. 
The stellar population forming after the accretion shows a 
lower helium abundance relative to the stars formed earlier 
from the pure ejecta of the super-AGBs. This may explain 
the presence of the very high-helium population observed in 
a few, very massive GCs. 

In Sect.|6]we present the results of a number of A^-body 
simulations exploring the subsequent cluster dynamical evo- 
lution. We show that the early cluster expansion driven by 
mass loss due to SN ejecta leads to a strong preferential loss 
of FG stars. Our simulations show that dynamical evolution 
can lead to values of the ratio of the number of SG to FG 
stars consistent with those suggested by observational stud- 
ies. We also show possible evolutionary routes leading to 
the loss of most of the FG stars, leaving an SG-dominated, 
cluster. 



2 THE PROGENITORS OF THE SECOND 
GENERATION: PREAMBLE 

2.1 Why AGB and super- AGB stars 

Although the study of the formation and dynamical evolu- 
tion of multiple population clusters is the main goal of this 
paper, one of the key assumptions of our investigation is 
that the SG progenitors are AGB stars, a scenario which has 
been criticized in a few recent studies aimed at producing a 
plausible chemical evolution history of GCs. In this section 
we briefly review a number of recent studies of the chem- 
ical modeling of SG stars and discuss why the AGB stars 
(along with the implications of the role of super-AGBs) are 
still to be considered a key ingredient in any theory for the 
formation of a SG population. 

In the last decade, attention has been focused mainly 
on the AGB scenario, based on the simple idea that the 
AGB winds are slow enough to be retained in the po- 
tential well of the cluster, and that these winds are pro- 
cesse d through the hot CNO cycle by "Hot Bottom Burni ng" 
( Ca meron fc Fowle j[l97ll : ISackmann fc Boothrovdlll992l') a t 
the bottom of the convective envelope ( Ventura et al.l200ll ). 
Many different groups have addressed the computation of 
suitable models for the chemical anomalies, and a num- 
ber of result s apparen tly do not su pport this scenario 
(|Fenner et al.i r2004: Ka rakas et al.|[20b6 : Choi & Yi 200^), 
while the model presented in iBekki et al.l (j2007f ) is in oifly 
partial agreement with the data. 

The model of SG formation by the contribu- 
tion of FRMS has been proposed much more r e cently 
(■Maeder fc Mevnetl I2OO6I : ]Prantzos fc Charbonnell I2OO6I : 



iDecressin et al. 2007al bh and, while it is not our goal to 
dismiss the FRMS scenario, we emphasize that many of its 
aspects are still to be fully explored. In particular, the main 
assumptions of this model are particularly stringent: all the 
stars must rotate at break-up speed, both to achieve a deep 
rotational mixing of CNO processed matter in the envelope 
and to inject, with small or zero ejection velocity, the CNO 
processed ejecta into the clu ster, where they mix w ith the 
intracluster pristine matter (IDecressin et al.l l2007ah . Con- 
versely, the high velocity winds during the helium-core burn- 
ing phase and the supernova ejecta must leave the cluster 
without any interaction with the remnant gas. The hydrody- 
namical and dynamical aspects of the FRMS scenario have 
not been studied, and a recent VLT FLAMES survey of mas- 
sive stars, studying the nitrogen abundance and the rota- 
tional velocity of 135 B-type stars in the LMC (|Hunter et al.l 
[200a), has found a number of unenriched fast rotators and 
cast some doubt on the efficiency of rotational mixing. 

As far as the AGB scenario is concerned, it is impor- 
tant to emphasize the following points. In discussing AGBs, 
we generally think of the low-mass, thermally pulsing AGBs 
that go through long series of thermal pulses and dredge- 
up episodes, so that they are generally rich in carbon and 
s-process elements. But these are not the AGB stars that 
can be progenitors of the SG in clusters. Some studies show 
that the sum of CNO elem ents is approximately constant 
a mong FG and SG star s (|Cohen. Brilev fc StetsonI |2002| : 
I Cohen fc Melendej [ioosh . This is possible only if we are 
dealing with the most massive AGBs, in which the episodes 
of third dredge-up are small in number, and the temperature 
at the base of the convective envelope is so large (> 8 x lO^K) 
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that the ON cycle is active, red ucing the O abundan ce as 
observed in the anomalous stars. ICarretta et alj (|2005t ) find 
that the CNO abundance is actually somewhat — although 
not much — larger in the SG st ars. The lat ter data indicate 
that the third dredge-up (e.g. Ilben fc Renzini 1983,) plays 
a (small) role, but the dredge-up episodes must be limited 
in number. Moreover, a few dredge-up episodes do not al- 
ter the s-process abundances, so that the lack of s-process 
enrichment in the SG is not an argument against the AGB 
polluters0 

A confusing issue when discussing the AGB evolution 
of low-metallicity stars is that the chemical yields based 
on AGB modelling by different groups differ enormously. 
The yields based on sets of models computed by Lattanzio's 
gr oup, widely used for f ollow- ups and recently published 
bv iKarakas fc Lattanziol l|2007l ). provide the correct abun- 
dance signatures of the SG only for the most massive AGBs 
(~ 7Mq ), while the yields of Ventura's grou p, and especially 
the most recent models bv lVentura fc D'Antonal (|2008al lbh. 
look compatible with the anomalies down to ~ SM© . The 
meiss-loss rate formulation, the nuclear reaction rates, and 
other details such as the treatment of core-overshooting in 
the phases previous to the AGB, are all relevant to th e 
determination of yields (e.g. IVentura" fc D'Antonal l2005bD . 
However, the most important difference among the physical 
inputs to these two sets is the treatment of superadiabatic 
convection, which affects the temperature stratification of 
models in the phase of "Hot Bottom Burning" and the ste l- 
lar luminosity and mass loss l| Ventura fc D'Antona|[2005al ). 
A "more efficient" convection model — like the FST model 
l|Canuto. Goldman, fc Mazzitellill 19961 ) adopted in Ventura's 
computations — destroys oxygen more efficiently, shortens 
the AGB lifetime, and reduces the number of thermal pulses 
and the effect of the third dredge-up for a wider range of ini- 
tial masses than in iKarakas fc Lattanziol (|2007| ). Thus, the 
AGB yields are very model-dependent, and the models of 
chemical evolution based on these yields are very model- 
dependent too, so that negative results can not necessarily 
be taken as falsification of the AGB progenitor scenario. 

The first complete model for the chemical anomalies by 
iFenner et al.l l|2004 ) shows that self-enrichment by AGBs is 
at variance with most of the observational anomalies, but 
this result could be easily predicted, as it is based on the 
IKarakas fc Lattanziol (120071 ) set. 

More recently, attention has been mostly devoted to 
explaining the presence of a population with extreme he- 
lium content, although, as we have noted, this population is 
only the tip of the iceberg of the SG, and it is present only 



^ A different, borderline, situation seems to be present in 
NGC 1851, w hose HR diagram shows a splitting of the sub- 
giant branch l lMilone et al.l I2OO8I ). that can be interpreted as 
due to two different stellar pop ulations differing in total CNO 
abundance (*Cassi si et al.l I2OO8I) . This is the only cluster in 
which a bimodal spread in the abundance of s-process elements 
has also been f ound, anticorrelated with the oxygen abundance 
llYone fc Grundahl.2008i ). More than others, this case might in- 
dicate that the AGB scenario is to be preferred to the FRMS 
scenario; with respect to most GCs, the SG of NGC 1851 may 
have been born from somewhat less massive AGBs, for which the 
effect of the third dredge-up is more clearly evident, both in the 
total CNO abundance and s-process elements increase. 



in the most massive clusters. As uj Cen sh ows the signa- 
tures of a complex chemical history (see e.g. I Romano et al.l 
l2007l . for a summary of the extensive observational studies 
of UJ Cen ) , and its formation took place most likely o ver an 
extended period (2-4 Gyr, e.g. iStanford et al.l [20061 ) ■ sev- 
eral attempts have been made to understand its very high 
helium popu lation by means of s tandard chemical evolu- 
tion models. I Romano et al.l l|2007h showed that the evolu- 
tion of UJ Cen can not be understood in terms of a "closed 
box" model, and that, also after relaxing this hypothesis, 
the helium-to-metal enrichment ISY / AZ ~ 70 of the blue 
MS stars of this system can not be attained, even with the 
very high helium yields of the winds of massive stars. Mod- 
els of UJ Cen taking into account inhomogeneous pollution 
by SNe la can reduce (by ~40%) but not eliminate the extra 
amount of helium need ed to explain the prese nce of the blue 
main sequence (bMS) (iMarcofini et al.ll2007l ). 

IKarakas et al.l (|2006l ) assumed that the helium yields 
come from the massive AGBs, and that the IMF is in- 
deed peaked at these masses, in order to maximize the ef- 
fect of the self-enrichment. In this case, they indeed could 
find a fraction of SG stars with Y up to ~ 0.35, but this 
is accompanied by a very large CNO enrichment, that is 
not observationally found. This result again depends on 
the models adopted to desc ribe the helium yield: only the 
IKarakas fc Lattanziol (|2007l ) models provided helium abun- 
dances as large as Y ^ 0.35 in the AGB ejecta, and this 
value was partially a chieved at the "second dredge-up" 
Ilben fc Renzinilll983l ). and partially by the third dredge- 
up during the thermal pulse phase. The third dredge-up of 
course also alters the total CNO abundance in the envelopes. 

Only recently has it become clear that a very high he- 
lium content (Y~0.38) can be achieved already after the sec- 
ond dredge- up in the env elopes of the most massive "super- 
AGB" stars l|Siessll2007al ) . The sum of the CNO abundances 
can remain close to the initial value during the following 
AGB thermal pulse phase, as the efficiency of third dredge- 
up is limited, becau se the helium luminosity during the ther- 
mal pulses is weak (ISiessll2007bl) . Unfortunately, full models 
through the super-AGB thermal pulse phase are not yet 
available, but it is clear that, in principle, it is possible 
to achieve high helium ab undances without e nhancing the 
CNO. FoUowing this idea, iPumo et all (|2008l ) suggest that 
the extreme helium-rich stars are formed directly from the 
ejecta of super-AGB stars. 

Finally, IChoi fc Yil | 20081 ). even adopting a model that 
maximizes the possible role of AGBs, tried and failed again 
to model the uj Cen /SY/ /S.Z ~ 70 for 30% of its population, 
concluding that "alternative processes are desperately called 
for." Again they do not include the most recent yields for 
super- AGBs in their computations, and use an almost closed 
box model, assuming that for the formation of the FG a 
fraction from 0.5 to 1 of the available gas is used. Thus 
the gas from which the SG is formed comes at most from 
half of the mass of the pristine gas plus the FG ejecta, or 
entirely from the FG ejecta. Under this hypothesis, for any 
reasonable shape of the IMF, the helium rich AGB ejecta can 
not reach 30% of the total cluster mass, as we will discuss 
again in Sect. 13. Il In fact, modelling must be able to predict 
at least a factor ten mass loss from the FG, to produce the 
desired number ratio of SG/FG stars of ~ 1. 
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2.2 The super- AGB phase 

In the exploratory models presented in this work, we do not 
consider in detail the peculiarity of the super- AGB evolu- 
tion, although we plan to include it in much more detail 
in subsequent work dealing explicitly with the chemistry of 
the ejecta. However, we discuss briefly here the super-AGB 
phase and its possible evolutionary outcomes, as it has con- 
sequences for our interpretation of these preliminary model 
results too. The minimum mass evolving into a core-collapse 
supernova is today referred to as Mmas , to distinguish it from 
the limit A/up, representing the minimum mass for carbon 
ignition, below which stars evolve into C-O white dwarfs 
(WDs). Stars in the mass range Afup< M <Afmas are called 
super- AGBs, and ignite carbon off-center in semi-degenerate 
conditions, but are not able to ignite hydrostatic Neon burn- 
ing in the resulting O-Ne core. Consequently, degeneracy 
increases i n the core, and these stars may undergo therma l 
pulses (e.g. Ilben et al.|[l997l : iRitossa et atll 19991 : ISiessll2006l ) 
and lose mass as "normal" (but quite massive and luminous) 
AGBs, the difference being the core composition (O-Ne ver- 
sus C-O) and core mass (^ 1.05 Mq ). 

As mentioned in Sect. 12.11 super-AGBs may play an 
important role in the evolution of proto-GCs. ISiesd (|2007al ) 
shows that their helium content in the envelope may reach 
the values Y=0.36 - 0.38, reasonably consistent with the 
highest helium content foun d in th e extreme populations of 
a few GCs, and lPumo et al.l l|2008l ) suggest that they are in 
fact the progenitors of the extreme He populations. 

The fate of super-AGBs depe nds on the competition 
between ma ss loss and core growth (|Poelarends et al]|2007l : 
ISiessI l2007ah : if mas s loss wins, t hey evolve into massive 
O-Ne white dwarfs (|Nomotol 11983 ) ; if the core reaches the 
Chandrasekhar mass, they evolve into electron capture su- 
pernovae (ecSNe), electrons being captured on the Ne nu- 
clei. Thus also stars below Mmas may explode as supernovae, 
but these events a re at least a factor t en less energetic than 
SNe II (^ 10^°erg |Dessart et al.lliooel ') and also a factor 10 
smaller than SNe la. In fact, it has been recently proposed 
that practically all the neutron stars (NSs) today present in 
GCs are born from ecSNe: due to their lower energy output, 
it is also probable that the newborn NS receives a propor- 
tionally sm aller natal kick, allo wing it to remain bound to 
the cluster jlvanova et akllioO^ . 

The presence of super-AGBs has an interesting conse- 
quence on the cluster evolution: the fate of super-AGBs in 
fact is either to end their life quietly as O-Ne WDs, or to 
explode as ecSN (either directly, or due to accretion induced 
collapse induced by mass transfer by a binary companion). 
In both cases, until O-Ne cores are formed, we do not have 
SN la explosions, that are due to the explosion of C-O WDs 
reaching the Chandrasekhar mass by accretion. These events 
can take place only after C-O WDs are formed in the clus- 
ter, when masses M <Mup evolve. Consequently, the epoch 
during which the super-AGBs evolve is probably the quietest 
period of the cluster's life, perturbed at most by ecSNe ex- 
plosions, either in single stars or by accretion induced col- 
lapse, and in any case not so energetic to alter either the gas 
evolution, as we will see in the hydrodynamical simulations, 
or its chemistry, as the whole core remains l ocked i nto the 
remnant neutron star. We refer the reader to ISiessI (|2007al ) 
for a careful study of A/mas as a function of the stellar metal- 



licity. His value for the typical metallicity of GCs showing 
multiple stellar populations is 10 Mq for models not in- 
cluding core overshooting. For models with moderate core 
overshooting this value is reduced by ^ 1 Afp) , thus to 9 Mq . 
The value of A/up from the models bv lVentura fc D'Antonal 
(200 8113) is - 6.3 Mq evolv ing at ~65Myr (see Table 1 in 
iVentma fc D'Antonal l|2008al )'). These latter models provide 
yields consistent with the GC chemical anomalies down to 
the evolution of 5 Mq , evolving at an age of ~ 10* yr. In 
this scheme, in summary, the SG may be formed by super- 
AGBs down to M> 6.3Mq , and by AGBs down to > 5Mq . 
In Sect. 14.41 we choose 40 Myr as the minimum age to start 
SN la explosions, in order to maximi ze their possible role. 
This lim it shifts to ^ 65Myr for the I Ventura fc D'Antonal 
(|2008al ) models, allowing a longer time for the formation of 
the SG. 



2.3 Post- AGB evolution 

We now consider whether the post- AGB evolution can play 
a role in the formation of the SG. In fact, after the loss 
of the envelope, the remnant hot cores evolve toward the 
white dwarf stage, reaching very high values of Tcff and thus 
becoming powerful sources of UV photons, acting as extra 
energy sources and potentially affecting the behaviour of the 
AGB gas. A study of the possible dynamical role played by 
extra energy sources is presented in Sect. 14.51 

Post- AGB evolution of the massive AGB core begins 
when the envelope mass is reduced below a critical value 
that does not allow stationary H-burning on the AGB, so 
that the residual envelope contracts (at the constant lumi- 
nosity dictated by the core mass itself) and the star becomes 
increasingly hotter. Planetary nebula excitation occurs at 
Tcff>30000K, when UV photons begins to be emitted. The 
critical envelope mass is a strong function of the core mass. 
In addition, the stellar luminosity is also i ncreasing with th e 
core mass. Based on the computations bv lPaczvhskil l|l97lh . 
Ilbeni l|l985h showed that the "fading time" that is required 
for the luminosity to fall by a factor 10 after the star has 
reached Toff=30000K is proportional to M^Jrl, so that, while 
it is >10000yr at the core masses of typical planetary nebula 
nuclei (O.GA/q), it is only a few hundred years for the core 
masses 1.0-1.3 Mq that are remnants of the possible pro- 
genitors of the SG we are considering (see also ■ Van Winckell 
l2003l . for a review). These objects, being so short-lived, are 
not a significant extra energy source capable of affecting the 
gas dynamics during the evolutionary phase we are consid- 
ering. 



3 INPUTS TO THE GAS DYNAMICAL 
MODELS 

3.1 The framework 

The model-building is as follows. The present GC is as- 
sumed to be composed of two stellar components, namely 
the FG and the SG stars. Its mass can thus be written as 

Afnow = AfpG,now + AfsG.now, where AfpG,now and AfsG.now 

represent the amount of mass today in long-lived FG and 
SG stars, respectively (recall, as mentioned in Sect. [TJ that 
we define long-lived stars as those having initial masses in 
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the range 0.1 < M/Mq < 0.8). The value ol x (in the range 
^ a; ^ 1) indicates the relative abundance of the two pop- 
ulations: Mfg .now - — xM-aaw and Msg ,now = (1 - x)Mnow. 

The FG is assumed to be already in place with an initial 
mass Mfg following a continous King radial profile 



1 



(1) 



up to the truncation radius rt. These stars have masses in 
the range 0.1 < M < 100 Mp distrib uted according to a 
Kroupa IMF $(m) (|Kroupa et al.ll 19931 ). The mass fraction 
5 of long-lived FG stars is given by 

/.0.8 

5 — m${m)dm. (2) 

J 0.1 

For our assumed IMF, we find 5 — 0.5. 

The evolutiona ry times are taken from 

iPadovani Matteuccil (|l993h . We will consider as pro- 
genitors of the SG the stars of masses below the limit for 
core-collapse supernova (SN II) explosion, here assumed 
to be M = 8 M0 . We discussed in Sect. [2]2] that, for the 
metallicities of GCs, this limit should probably be shifted 
to ~ 9 Mq for models including core overshooting, but the 
precise choice does not make an important diff'erence to the 
global picture. 

We start our calculations at f = 28 Myr, when all 
the SNe II have exploded, clearing the GC of all the pris- 
tine gas (we will reconsider the possible role of pristine gas 
in Sect. [Sjl. At this stage, Mfg is reduced by an amount 
AMfg = eMfg; for a Kroupa IMF e = 0.08, but it can be 
much larger for flatter IMFs (see Sect. I4.3|l . The SG stars 
form from the ejecta of the AGB stars of the FG during 
a period of time Atf, which must be long enough to allow 
the formation of a substantial SG population, but not so 
long as to include evolution of stars of masses M^C 4 — 5Mq, 
which do not provide yields compatible with the observed 
chemical anomalies. The value adopted in our simulations, 
Atf = 100 Myr, satisfy this constraint and only stars with 
initial masses M > 4M0 contribute to the SG formation. 

The amount of gas returned by the FG in this period 
is A Mg = PMfg, where, for Kroupa IMF (3 = 0.05. Fol- 
lowing JPrantzos fc Charbonnell (|2006l ) and iDecressin et al.l 
l|2007al ). we assume that SG stars have masses within the 
range 0.1-0.8 Mq, and thus are all still present today. In the 
working hypothesis that all the gas forms stars, i.e. that the 
star formation efficiency is r; = 1, the SG mass is therefore 
MsG.now = AMg = /3Mfg. 

We stress that our assumption of a SG IMF skewed to- 
ward low masses is made to minimize the FG mass: if stars 
with M > 0.8 M0 are assumed to also be formed, the re- 
quired ejecta mass, and thus Mfg, would be larger. The 
assumption of a flatter IMF, with its larger mass return, 
does not lead to a significant lowering of A/fg. In fact, in 
this case, the number of long-lived FG stars would decrease, 
and thus a large value of Mfg is still needed to obtain a sig- 
nificant quantity of them (see Sect. 14.3) . It is worth noticing 
that while SG stars more massive than 0.8 Mq can be in- 
cluded in our model, stars with masses larger than 8 Mq 
must be excluded; these stars explode as SNe II and would 
sweep out the gas returned by the FG stars, preventing the 
formation of a substantial stellar SG. In any case, we point 
out that formation of (only) low mass stars has been invoked 



also in the cooling flow occ urring in elliptical galaxies (e.g. 
iMathews fc Brighentilll999l ) which are a scaled version of 
the cooling flows taking place in our GC models (see Sect. 

From all the above, it follows that Mfg = /Mnow, 
where / = (1 — x)//3. For a Kroupa IMF, and assuming 
X = 0.5 as a rule of thumb, we obtain Mfg = lOMnow. 
Thus, a GC with a current mass Mnow originates from a FG 
ten times larger, or even larger if the gas is not completely 
turned into stars (or if the IMF of the SG is not limited 
to 0.8 Mq). This huge extra amount of mass must be lost 
in the successive evolution of the cluster. We postpone the 
discussion of this point to Sect.[Sl and describe in the next 
sections the dynamics of the gas returned by the FG and 
the formation of the SG stars. 



3.2 Physical inputs 

In order to follow the gas evolution in the GC, we 
adopt the ID ver si on o f the hydrocode described in 
iBedogni &: D'Ercold (|l986l ') and solve numerically the fol- 
lowing equations: 



dp 
dt 



+ V ■ (pv) = ap. 



dm „ , , , m 

— + V ■ (m (g) v) = -(7 - 1)V£ -sp-u — 

Ot tgf 



(3) 
(4) 



^ V ■ (_Bv) = -(7- 1)SV • V- - i + S (5) 
at tsf 

In the above equations, v, p, m and E represent, respec- 
tively, the gas velocity and the densities of the gas mass, 
momentum and internal energy. The gravitational acceler- 
ation g is due to both FG and SG stars, as well as to the 
gas. The ratio of the specific heats is 7 = 5/3. Given the as- 
sumptions described above, the mass loss occurs only from 
FG stars at the specific rate a = Mfg /Mfg, and p« thus 
represents the FG stel lar density. The ma ss loss rate Mfg 
is computed following ICiotti et al.l (|l99ll ). adapted to the 
IMF chosen for the present case. The star formation rate 
(SFR) is characterized by the timescale 4f ~ niax(tcooi, iff), 
where tcooi and tg are the cooling time and the free-fall time 
of the gas, respectively; it is usually tcooi <C iff and thus 
tsf oc iff (X p""'^. The arbitrary parameter u summarizes 
our uncertainties about star formation, but we find that our 
results are essentially independent of its value (the reason 
for this is discussed in the Appendix); as a consequence, all 
the models presented in this paper have 1/ = 1. We stress, 
however, that the star formation is allowed only in regions 
where the gas is converging. For this reason we set 



0, V ■ V ^ 

1, V-v<0 



(6) 



Radiative losses are given by L = n'^A(T), where n 
and T are the number density and the temperature of the 
ga s, respectively, and A(r) is the cooling curve as given 
bv lRosen fc BregmanI l| 19951 ). The source term S takes into 
account the heating of the gas due to the feedback of the 
stars: 



S = 0.5ap. (3a 



2q) 



(7) 



© 2008 RAS, MNRAS OOO.nifSTI 



Multiple generations in globular clusters 7 



The first three terms in the above equation are due to the 
thermalization of the kinetic energy of the AGB winds, 
where i; is the gas velocity, a is the one-dimensional ve- 
locity dispersion of the FG stars, and v„ = 2 x 10^ cm s~^ 
is th e wind velocity of the AGB stars (e.g. [Marshall et al.l 
I2OO4I ) . We also computed a set of models in which an ad hoc 
heating source q of stellar origin is added (see Sect. 14. S) . 

Assuming that the SNe la belonging to the FG start 
to explode with a delay of 80-100 Myr, we are allowed to 
neglect their effect in our simulations, and Eq. [7] represents 
the total energy source of the system. If, however, such a 
delay is shorter, the repercussions of the stellar explosions 
must be taken into account (see Sect. 14.4")) . 

We conclude this section by pointing out that in our 
hydrodynamical simulations we have assumed a static stel- 
lar potential. As discussed in Sect. [6] the cluster is expected 
to expand during its early evolution in response to the early 
mass loss due to SN ejecta. In a further study now underway, 
we will explore the role of the early expansion of the cluster 
driven by primordial gas loss and SN II ejecta on the for- 
mation and final structural properties of the SG system. We 
have carried out a number of preliminary simulations with 
initial conditions as in the standard case, for a FG system 
expanding homologously and in isolation; the results seem 
to indicate that the central gas and SG stellar densities and, 
in general, the conclusions concerning the SG properties for 
the standard case are not significantly affected by the cluster 
expansion. 

3.3 The energy budget 

The energy required per unit time to extract f rom the GC 
the g as returned at rate a by the FG is (see ICiotti et al.l 

r 2 

Lg = -a{t) I inr pFG{r)4'FGir)dr, (8) 
Jo 

where pfg('') and (/f>FG(r) are the radial profiles of the FG 
stellar density and gravitational potential, respectively. The 
energy input linked to such mass return and due only to the 
stellar velocity dispersion is 

i<T = |a(t)^ 4nr'^PFG{r)a^{r)dr. (9) 

By definition, these two equations can be written as 

Lg^2a(t)\U\, (10) 

L„ = a{t)K, (11) 

where U and K are the binding energy of the GC and the 
kinetic energy of the stars, respectively. If the FG is in virial 
equilibrium, then Lg — ALa. This is a general result inde- 
pendent of the amount of mass and of the shape of the FG 
potential well. The first three terms in Eq. [7]are of the same 
order of magnitude; therefore in principle an extra source 
Q ^ I/O- would be enough to revert the cooling flow into 
a wind, thus halting the SF. However, because of the large 
radiative losses occurring in a cooling flow, a wind can start 
only for sources much larger than Lo- (see Sect. 14.5) . It must 
also be pointed out that, for sufficiently small Mfg, the en- 
ergy of the AGB winds may vent from the GC at least part 
of the gas returned by the stellar FG (see Sect. I4.2 |l. 



4 RESULTS 

4.1 The standard model 

For our standard model we assume a FG stellar mass Mfq = 
10^ M0 in order to describe a GC with a current mass 
AInow ~ 10^ Mq , as discussed in Sect. 13. II This initial stellar 
population follows a King radial profile with the parameters 
below (see Eq. [T|: po = 10^ pc~^, = 6.3 pc and 

rt = 200 pc. The truncation radius of the model coincides 
with its tidal radius at a distance of 4 kpc from the Galactic 
center. 

Fig. [T] shows the radial profile of several quantities at 
t — 100 Myr, when the mass of the SG stars has grown to 
Msg = 4.45 x 10^ Mq . The upper panel shows that the fluid 
velocity is negative over the entire cluster; radiative losses 
prevail over the heating source, particularly in the central 
region where the amount of gas shed by the FG is larger 
and causes significant radiative cooling, which strongly in- 
creases with the gas density (see Sect. [5^ . Consequently, 
the gas pressure falls and gravity succeeds in driving the 
gas inward, establishing a cooling fiow from the start. This 
gas accumulates in the center, as shown by the thick solid 
line in the lower panel of Fig. [1] The SG stars forming from 
this gas are thus more concentrated than the FG stars; their 
distribution is indicated by the dashed line, terminating at 
radius r ~ 100 pc, half of the assumed rt. The lower panel 
also shows the temperature profile increasing with radius, 
as expected in a cooling flow. 

Since tst oc p~^'^ (see Sect. 13.2) . the SFR is proportional 
to p^'^ (see Eq. [3]), and the SG stars form preferentially in 
the central region, where they are highly concentrated. The 
dashed line in the lower panel of Fig. [T] indicates the radial 
distribution of the SG stars as they form from the gas. 

The thick lines in Fig. [2]show the gas and the SG mass 
evolution, as well as the time evolution of the SFR. Although 
the star formation process is rather effective, we point out 
that AfsG = 5 X 10^ M© (corresponding to the case x = 0.5) 
is attained only a.t t = 120 Myr. Actually, the flnal value of 
the SG mass depends on the time at which a suitable number 
of SNe la belonging to the FG start to explode, removing 
gas from the cluster and inhibiting further star formation. 
We will discuss the role of SN la explosions further in Sect. 



4.2 Model with a lower initial FG mass 

We have run a model similar to the standard one, but with 
a stellar FG mass A/fg = lO'' Mq. In this case, the initial 
parameters are po ~ 1.08 x 10^ M© pc~^, ~ 2.85 pc and 
rt = 90 pc (also in this case the truncation radius corre- 
sponds to the tidal radius at a distance of 4 kpc from the 
Galactic center). In principle, since the gravitational poten- 
tial, mass return rate, and La all depend linearly on A/fg, 
the gas behaviour in this case should simply be a scaled 
version of the cooling flow developed in the standard case, 
with lower velocities, temperatures and densities. However, 
the source term S (see Eq. (Tjl depends also on v^, the ve- 
locity of the AGB winds, which does not scale with Mfg- 
Moreover, radiative losses depend quadratically on the gas 
density (see Sect. 13. 2p . For these reasons, the self-similar 
scaling of the models with Mfg breaks when Mfg is low 
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Figure 1. Radial profile of several quantities at t = 100 Myr 
for the standard model . Upper panel: gas velocity (solid line) 
in units of 10® cm s~^ and pressure (dashed line) in units of k 
erg cm~^, where k is the Boltzmann's constant. Lower panel; gas 
density (thick solid line), gas temperature (thin solid line), SG 
stellar profile produced by the simulation (dashed line). 
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Figure 2. Evolution of the standard model. Upper panel: SFR 
evolution. Lower panel: temporal profile of the amount of gas 
(solid line) and SG stars (dashed line). The thick lines refer to 
an evolution in absence of SNe la. In both panels the thin lines 
indicate the effect of SN explosions starting at t = 40 Myr and 
occurring every 5 X 10* yr. 
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Figure 3. As Fig.[T] but for Afpc = lO'' M© 
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Figure 4. As Fig. [2] but for AfpG = 10*' Mq . 



enough to generate a stellar velocity dispersion a a factor 
of a few lower than and a mass return from the FG 
stars such that the gas density is not so large to provide 
effective radiative cooling. In the present model the source 
term prevails on the radiative losses only at the outskirts 
of the cluster; as a consequence, a partial wind forms, with 
the stagnation point located at r = 23 pc (see Fig. [Sjl . For 
this reason the SG stellar population is relatively more con- 
centrated than in the standard case, extending only out to 
r = O.lSn. However, the comparison between Fig. |2] and 
Fig. |4] (thick lines) shows that the global behaviour of gas 
and SG stars are similar in the two models. 
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Figure 5. As Fig. (2] but for a flatter IMF (see text). 

4.3 Model with a flatter IMF 

FoUowine iD'Antona fc Caloil ()2004h . one can assume a flat- 
ter IMF for the FG in order to get a larger number of giant 
stars at early epochs. In this way, the larger amount of gas 
returned by such stars would facilitate the formation of the 
SG. To explore this scenario, we assume that for M < 5 
M0 the IMF follows a power law with an exponent s — ~1, 
while more massive stars follow a Salpeter distribution with 
s = —2.35. The flatter IMF leads to a larger amount of Mg 
and Msg at any time compared to the standard model, as 
illustrated in Fig. [S] (thick lines). 

Repeating for this model the simple analytical calcula- 
tion of Sect. 13.11 we obtain f3 = 0.2, 5 = 0.05 and e — 0.5. 
Assuming x — 0.5, i.e. that half of the MS stars with 
0.1 < M/Mq < 0.8 currently in the GC belong to the FG, 
we have 0.05Mfg = 0.5Mnow and thus A/fq = lOAfnow. 
As AAfg = 0.2Mfg = 2Afnow, the star formation efficiency 
must be = 0.25. In conclusion, also in this case, if one 
requires the current mass of long-lived SG and FG stars 
to be similar, the total FG initial mass must be ten times 
larger than the current total mass of long-lived stars. How- 
ever, in this case the cluster must evolve without losing mass 
because the initial number of long-lived FG stars is much 
smaller. We also point out an important difference con- 
cerning the relative abundance of stellar remnants (white 
dwarfs, neutron stars and black holes) and long-lived MS 
stars (0.1 < M/Mq < 0.8). Speciflcally, in both cases at 
t ~ 13.5 Gyr, Mrom ~ 0.15AfFG, but for systems with a 
flatter IMF the total mass of long-lived stars is smaller and 
the ratio of the total mass in stellar remnants to the total 
mass in long-lived stars is much larger (by a factor ^ 10^) 
for systems with a flatter IMF. 

4.4 Models with SN explosions 

In order to understand the effect of SNe la in our standard 
model, we simulated the expansion of SN remnants (SNRs) 
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Figure 6. Velocity, temperature and density of the gas in the 
GC at At = 10'* yr (dashed line). At = 3 X 10* yr (long-dashed 
line) and At = 5 X 10* yr (dot-dashed line) after a SN la explo- 
sion occurred in the standard model at t = 40 Myr. The solid 
lines show the radial profiles at t = 40 Myr just before the SN 
explosion. After At = 5 X 10* yr another SN la is expected to 
explode (see text). 



through the GC. Given the spherical symmetry of the model, 
we assume that the explosions occur in the center; such an 
assumption is reasonable, as in the central region there is a 
larger density of SN la progenitors belonging to the FG. We 
also assume that the flrst explosion takes place at the time 
tpk at which the SN la rate quickly reaches its maximum. 
This time is not unambiguously determined, and its value 
for a SF of shor t duration is estimated in the range 40 < 
tpk < 300 Myr l|Matteucci fc Recchil [20011 : iMannucci et all 
'200Ci). In order to emphasize the effects of the SNe la, in the 
present models we adopt tpk = 40 Myr, f ollowing the SN la 
rate fllustrated bv lMarcolini et all l|2006l) . 

We flrst consider the explosion of a single SN la taking 
as initial values of the variables those obtained in the stan- 
dard case at t = 40 Myr. The occurrence of the supernova 
is mimicked by adding an extra amount _Esn = 10^^ erg of 
thermal energy in the flrst three mesh points of the grid. 

We point out that the binding energy of the gas at the 
moment of the explosion is only 44% of that of the SN, so in 
principle the GC could be completely deprived of gas by the 
stellar explosion. However, given the large density of the gas 
itself, the radiative losses are rather effective, and in fact no 
gas is lost by the GC. 

Fig. [6] shows the SNR expansion at three different 
times. As expected in a Sedov-like solution, a "hole" is 
present behind the shock front. Contrary to a genuine Se- 
dov solution, however, this hole is partially replenished by 
the FG high mass return rate that at this stage is nearly 
8.9 X 10~^ M0 yr~^. The gas replenishment favours the ra- 
diative losses, thus weakening the SNR shock and slowing 
the remnant expansion below the escape velocity. 

Although the effect of a single SN la explosion may 
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therefore appear as a minor perturbation of the coohng flow, 
the cumulative effect of numerous SNe la drastically alters 
the gas dynamics. At t = 40 Myr the SN la rate is of the 
order of 2 X 10"^ yr"^ for Mfg = 10^ M© (|Marcolini et all 
[2OO6), i.e. one SN explosion every 5 x 10* yr (the maximum 
time span considered in Fig. [Sjl . More recent SNRs overlap 
older ones before these latter merge with the medium, and 
the sum of their actions has substantial consequences on the 
GC evolution. As shown in Fig. (2] the SFR quickly declines 
and the SG mass stops to grow attaining the value Msg = 
8.4 X 10* Mq. This result is obtained assuming a constant 
SN la rate up to 100 Myr; this assumption is largely justified 
as the actual rate suffers only a negligible reduction during 
this interval of time. 

We point out that such a dramatic reduction of the SFR 
is due not only to a decrease in the gas content (the SFR 
being proportional to p^'^), but also to the condition for star 
formation we have adopted (see Eq. [SJ. In fact, the gas is 
set into expansion over much of the GC volume by the SN 
explosions, and dynamical conditions are not favourable to 
further star formation. 

In the end, as Msg is only one-fifth of the value attained 
in the standard case, in order to obtain a cluster with a 
current total mass Mnow ~ 1.7 x 10^ Mq with an equal 
amount of FG and SG stellar components, an initial FG 
population with Myg = 10^ M© must lose ~ 99% of its 
mass. 

It is interesting to note how differences in Mfg and/or 
the slope of the IMF influence the effectiveness of SNe in 
halting the formation of the SG stars. Fig. [5] shows that, 
contrary to the standard case, the flatter IMF essentially 
neutralizes the effect of the SNe la, as Msg is only very 
slightly perturbed by the stellar explosions. In fact, when the 
SNe la start to explode the amount of gas is Mg — 7.6 x 10* 
Mq, a factor of 2.1 larger than in the standard 
consequence the radiative losses are roughly 2.1^ = 4.2 times 



greater, and suppress rather effectively any disruptive effect 
of the SNe. 

Models with A/fg ~ 10^ M© show similar behaviour 
(see Fig. H] and Fig. [7J , independent of the IMF adopted. 
In this case the gas has time to recover, at least in part, its 
inward convergent flow between successive SN explosions, 
because the SN la rate is ten times lower than in the models 
with Mfg = 10^ M© . Despite their lower rate, however, the 
stellar explosions are equally effective in reducing the gas in 
the cluster because of its lower binding energy in this case. 

In principle, models in which the SF proceeds also dur- 
ing the SN activity would show today some stars polluted 
by the SN la ejecta. For example, let us examine the model 
with A/fg = 10^ M0 and flatter IMF, which is the most 
favourable case to lock SN ejecta in SG stars, as illustrated 
in Fig.O The mass return rate at 40 Myr is 3.4 x 10"^ M© 
yr~^, while the mean production rate of metals by SNe la 
is 2.8 X 10~^ M0 yr~^; as a consequence, the gas accu- 
mulating within the GC has a metallicity increment due to 
the SNe la AZ — 8.2 x 10~*, which is comparable to the 
typical metallicity of GCs showing self-enrichment. On the 
other hand, star to star iron variatio ns in GCs do not ex- 
ceed ~ 0.04 dex (jCratton et al.|[2004l 'l. so that retention of 
SN la ejecta in SGs is r uled out by observati ons, with the 
exception of ui Cen (e.g. iMarcolini et al.ll2007l ). 



4.5 Models with an extra energy source 

In this section we explore the effect of an extra energy source 
due, for example, to X-ray binaries or a number of luminous 
hot planetary nebulae nuclei. In our model, this source is 
distributed as the FG stellar density and its total amount is 
equal to 

Q= 47rr^pFG(r)gdr (12) 
Jo 
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Figure 9. Q^r as a function of the cluster concentration c for 
models with the same mass, IMF and truncation radius as the 
standard model. 
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We have carried out a number of exploratory simula- 
tions and found that, in general, a critical value Qcr exists, 
separating models hosting a cooling flow from models host- 
ing a wind. The transition between the two kinds of solutions 
happens rather abruptly as Q increases from values Q < 
to values Q > Qcr- Such a sudden transition is due to the 
strong nonlinearity of the radiative cooling term, which de- 
pends on the square of the gas density. 

Fig. |S]is analogous to Fig. [5] for a model similar to the 
standard one, but with Q — Qa = 6 x 10''^ erg s~^. After 
an initial transient phase during which some SF occurs, a 
steady wind is established which prevents any further sub- 
stantial growth of Msg- This case is borderline between a 
model with Q = 4 x 10^^ erg s"'^, in which an almost normal 
star formation occurs, and one with Q — 1.3 x 10"^** erg s^^, 
where the immediate onset of a fast wind precludes any SG 
formation. For this set of models, the role of SNe la is neg- 
ligible, as they start to explode when the SG star formation 
has already been halted. We stress that, for the standard 
model, we have = 2 x 10^'' erg s~^ at f = 28 Myr, much 
lower than Qcr, as expected owing to the radiative losses 
(see Sect.[33J. 

We mention that Qa = 2 x 10^^ erg s~^ for the model 



with Mfg = 10" 
IMF, we have Qcr 



M 





6 x 10' 



while, for the models with flatter 



erg s"^ and Qcr = 2 x 10"" erg 
s " for MpG = 10^ Mq and A'/pg = lO'^ Mq , respectively. 
As expected, Qcr is larger for flatter IMFs, because of the 
larger mass return rate which increases the radiative losses. 



4.6 Models with different initial concentrations 

In our standard model the FG is distributed following a 
King proflle with concentration c = log(rt/rc) = 1.5 (see 
Sect. We analyzed also cases with the same values of 
MpG and rt but different concentrations, namely c — 0.75, 
c = 1, c = 2 and c = 2.25. For these latter models, at t = 28 
Myr one has = 0.86 x 10^^ erg s~\ = 1.1 x 10^^ erg 
s"\ La = 3.5 X lO^'' erg s~^ and L,, = 5 x 10^^ erg s~\ 
respectively. 

Within the range of values of c we have considered, the 
results are remarkably similar to those of the standard case, 
and in fact Qcr varies only by a factor of ~ 2, as shown in 



logdv 



Figure 10. Radial profile of the SG stars at t = 100 Myr for 
models with MpQ = lO'^ Mq , c = 1.5 and rt = 40 pc (long- 
dashed line), rt = 100 pc (dotted line), rt = 200 pc (solid line, 
standard model), rt = 400 pc (short-dashed line). The masses are 
normalized to their respective maximum values: Msg = 4.8 X 10^ 
Mq (long-dashed line), Msg = 4.8 X 10^ Mq (dotted line), 
Msg = 4.45 X 10^ Mq (solid line, standard model), Msg = 
3.52 X 10^ Mq (short-dashed line). 



Fig. [9] There is a trend of increasing Qcr with c due to the 
fact that the stellar density gets larger in the GC core; as a 
consequence, the gas density increases because of the larger 
mass return, and the radiative losses are stronger. A partial 
wind establishes in the outer regions, while a cooling flow 
persists in the core. Reverting such a 'mini-inflow' becomes 
increasingly more difficult for larger values of c. We point out 
that, although La increases faster than the critical source, 
it remains always much lower than Qcr. 



4.7 Models with different initial truncation radii 

In order to understand the inffuence of the assumed size of 
the FG population on the relative concentration of the SG 
stars, we ran models with the same mass Mpg and concen- 
tration c of the standard model, but with different values of 
the truncation radius (vt = 40, 100, 400 pc). The compari- 
son among the models is summarized in Fig. [10] and shows 
only a minor dependence of the relative concentration of the 
two populations on the initial size of the cluster. In all cases 
the SG population is strongly concentrated in the cluster 
innermost regions, with 80-90 percent of the total SG mass 
enclosed within r/rt ~ 10""' — 10~^. 



4.8 Consequences of the end of the cooling flow 

The energy input either by recurrent SN la explosions or 
by distributed stationary sources with total power > Qcr 
drastically limits the epoch of gas accretion and SG star 
formation in the cluster core. Clearly, there are a number of 
uncertain parameters in our understanding of SN la rates, 
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and also their starting age and actual rate at early epochs 
may be smaller than we suggest (see Sect. 12.2)1 . 

Our simulations show that SN la explosions eflFectively 
clean clusters with initial masses in the range 10'^ - lO'' Mq 
from AGB ejecta , and limit the SG formation time. We 
point out that this is no longer the case for more massive ob- 
jects like the ultra-compact dwarf galaxies, and possibly the 
progenitor of lj Cen , which are embedded in massive dark- 
matter haloes, and which proceed to a mo re complex (but 
also more s tandard) chemical evolution (see lMarcolini et al.l 
I2OO6I . I2OO7I ). 

Concerning the stationary sources, Qcr ~ 10^* erg s~^ is 
of the order of the Eddington accretion luminosity for NSs, 
so a number of sub-Eddington accreting NSs distributed 
in the cluster could succeed in stopping the SG formation 
epoch. We know that several hundreds of NSs, revealed by 
their millisecond p ulsar radio emissi on, are present in to- 
day's clusters (e.g. Heinke et al] |200d). Their radio emission 
(relevant for those NSs that are already spun up to millisec- 
ond periods at the early epoch we are considering) is in the 
plausible range 10^* — lO'^^erg s~^, so it has no global effect 
on the cooling flow. Only accreting NSs are good candidates 
as stationary source. Although only a few stationary X-ray 
sources are currently active in all the galactic GCs, the situ- 
ation might have been different in the early stages of cluster 
evolution and this issue will need further investigation. As 
discussed in Sect. 12.21 NSs are just born from the highest 
mass super- AGBs, when the super- AGB winds are collect- 
ing in the core. 

From the chemical point of view, the fact that the cool- 
ing flow is halted is important since it prevents less massive 
AGB stars from contributing to the formation of the SG. 
Although the maximum allowed age (and minimum allowed 
mass) strongly depends on the efficiency of convection, third 
dredge-up and mass loss in AGB stars, if the cooling flow is 
not halted and star formation prolongs beyond ~ 10* yr (the 
upper limit adopted in the standard hydro models of Sect. 
I4.1|l . in all the current stellar models the third dredge-up is 
important, and the yields are not compatible with the ob- 
served chemical anomalies; specifically, the total CNO abun- 
dances and the sodium abundance will be much larger than 
observed, and the oxygen abundance will not be reduced. 



5 PRISTINE GAS 
5.1 Motivations 

If the SG formation time is too short, it will not meet the ob- 
servational requirements, as not enough matter will be con- 
tained in the ejecta. One possible solution is the following: 
while the SN II events have wiped out the pristine matter 
from the cluster center, during the AGB epoch it is pos- 
sible that the pristine gas in the outskirts participates in 
the cooling flow. In the central regions it will mix with the 
AGB gas and contributes to star formation. In this way, the 
span of SF time required to achieve a massive SG is shorter. 
Further, the "dilution" of the AGB ejecta by pristine gas 
may fully explain, in particular, th e milder chemical anom a- 
lies (iPrantzos fc C harbonncl 200^; iDecressin et al.ll2007al lbl: 
IVentura fc D'Antona.2008a .b). 

In all the models discussed in the previous sections, it 
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Figure 11. As Fig. [2] but for the infall model. The long-dashed 
line represents the amount of pristine gas. 

has been assumed that the FG SNe II deprive the GC of all 
its pristine gas. In this section we drop this assumption and 
study a model in which part of the pristine gas contributes 
to the SG formation. We assume that initial asymmetries 
in the distribution of gas allow to vent out the SN II ejecta 
along preferential directions, creating an "hourglass" cavity, 
and leaving part of the pristine gas in a substa ntially unper- 
turbe d torus at the outskirts of the cluster (e.g. lRecchi et al.l 
[200I). 

5.2 Model with pristine gas infall 

Since we use a 1-D code for our hydrodynamical simulations, 
we can not properly explore the case of SN II ejecta blown 
away through an axisymmetric chimney, with part of the 
pristine gas confined around it. We approximate this situa- 
tion assuming an initial gas distribution within the GC with 
a central hole of radius r — 150 pc and with a density radial 
profile p(r) oc beyond it. The total gas mass within the 
numerical grid is M = 2.6 x 10^ M©. These values are rather 
arbitrary, but nevertheless the model is useful to understand 
the consequences of the primordial gas infall on the SG for- 
mation process. All the other details of the model are the 
same as in the standard model. 

As apparent in Fig. 1111 after 100 Myr the total mass 
of SG stars formed, Msg = 6.2 x 10^ Mq, is larger than 
that of the SG population formed in the standard model. 
This larger amount is created very quickly because part of 
the pristine gas is in place from the beginning and available 
for star formation. As this gas moves toward the center, and 
more gas enters the GC from outside, the SFR increases. 
After ~ 50 Myr no more pristine gas falls within the GC, 
and that already present has been locked within the SG stars 
formed till then (see the dashed line in the bottom panel of 
Fig. Ill|l . As a consequence, the global evolution of the GC 
recovers its "standard" behaviour (see Fig. [2]). 

Additional interesting details of this model are illus- 
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Figure 12. Radial profile of several quantities of the infall model 
at two different times: total gas density (thick solid line), pristine 
gas density (long-dashed line), AGB ejecta density (thin solid 
line), SG stellar density (dashed line). 



trated in Fig. 1121 where the radial profiles of several quan- 
tities are displayed at two different times, namely t = 40 
Myr and t = 45 Myr. After 40 Myr the pristine gas has 
already largely moved toward the GC center. However, the 
GC central volume included within r 5C 10 pc still hosts only 
gas ejected by AGB stars, and the SG stars formed in this 
volume till this time have thus the same chemical character- 
istics of these polluters. Successively, up to t ~ 50 Myr (see 
Fig. Ilip , the pristine gas occupies all the GC volume and is 
well mixed with the gas returned by the FG stars. The stars 
forming in this period sufi'er a "diluted" pollution. Finally, 
for t > 50 Myr, the FG stellar ejecta are again the only 
source of gas for the SG star formation and the chemical 
composition of the new SG stars formed after this time is 
that of the FG progenitors' ejecta. 

At t = 40 Myr the SG stars formed at r < 10 pc have a 
mass ~ 8.5 x lO^M© , and they are polluted by AGB stars 
with masses M ^ 7 M© which are strong He polluters. For 
40 < t < 50 Myr the chemical contamination of the SG stars 
is diluted by the pristine gas all over the GC, and is due to 
FG stars with masses of 6-7 Mq . 



5.3 Implications of the pristine gas infall for the 
very helium-rich population 

Although the study of models including the infall of pristine 
gas will need to be extended and further explored, a number 
of interesting results have already emerged from the prelim- 
inary simulation discussed in the previous section. 

We have shown that there is an epoch, between the 
end of SN II phase and the time when the pristine gas be- 
gins to mix with the super-AGB gas ejecta, during which, 
if the cluster is massive enough, the winds from the super- 
AGB stars are the only contributors to the cooling flow and 
the SG formation. We suggest that this is the phase dur- 



Figure 13. Histogram of the helium abundance of the SG stars 
formed in the model with the pristine gas infall (solid line). The 
dashed line shows the histogram of Y for NGC 2808 inferred by 
D'Antona and Caloi (2008) from the observations of the HB by 
Bedin et al. (2004) and from the MS distribution by Piotto et al. 
(2007). Each histogram is normalized to its maximum value. 



ing which the most extreme very helium-rich SG population 
can be formed. The super-AGBs can indeed have Y~0.38 
in their ejecta, and their gas will not be diluted. Thus a 
very homogeneous population with very high helium will be 
formed. 

This early phase of very helium-rich SG star formation 
is followed by an epoch during which the pristine gas reaches 
the cluster inner regions and the SG stars form from gas in 
which the helium of the AGB ejecta, collected in the inner 
regions by the cooling flow, is diluted with the primordial 
helium. If the cooling flow goes on after the pristine gas is 
exhausted, stars will again be formed with undiluted ejecta, 
but the evolving AGB mass and its Y will be smaller (e.g. 
IPumo et al.ll200a ). giving rise to a "helium gap" between 
the first very high helium stars and the other SG stars. 

The scenario outlined here might lead to the formation 
of the blue main sequences in lu Cen and in NGC 2808, 
and of the extremely blue HB stars in o ther clusters (for a 
parametric approach to this problem, see iBekki et al.ll2007l ). 

We implemented in the code the helium abundance of 
the eje cta as a function o f the initial stellar mass, follow- 
ing the IPumo et al. ll2008h comp i lation, based on the AGB 

r^T 1^ til ci Viir I X/Qn + i 1 T-o 



models by IVentura & D'Antona' (^20083) and on the super- 
AGB models by Siess (2007a). Fig. 1131 shows the predicted 
SG helium distribution for the infall model. 

We see that a peak in the mass distribution is obtained 
at Y~ 0.37, corresponding to the star formation in the cool- 
ing flow preceding the infall of pristine matter, and another, 
well distinct broader peak at medium-He (Y~ 0.31 — 0.34) is 
obtained from star formation when the wind matter is mixed 
with the pristine matter. There is a negligible tail of stars 
formed at Y< 0.25, so that this model implies that the 50% 
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stars having normal helium in NGC 2808 must indeed be a 
remnant of the FG population. In Fig. [13] we also show the 

Y distribution derived for th e SG of the cluster NGC 2808 
bv lD'Antona fc Caloil (|2008l ). 

Although we have not made any attempt to find the 
initial parameters leading to the best fit of the observational 

Y distribution, the theoretical and the observed distribution 
are, in general, remarkably similar. 

Since the helium yields from super- AGBs and AGBs are 
still uncertain, the amount of pristine matter involved, and 
its mixing with the FG matter is yet to be fully explored, we 
believe our results show that this scenario is a very promising 
step towards the understanding of the origin of the super-He 
rich population. 



6 STELLAR DYNAMICAL EVOLUTION 

One of the key requirements of the standard scenario pre- 
sented in the previous sections is the large initial mass of 
the FG population. This assumption, as we have discussed, 
is necessary for the formation of a significant SG population 
from the ejecta of intermediate-mass stars. The next issue to 
address in order to test the viability of the proposed scenario 
concerns the possible mechanisms driving the evolution of 
the ratio of the number of SG to FG stars. 

Numerous s tudies of the dynam ical evolution of star 
clusters (see e.g. iHeggie fc Hutll2003l . and references therein 
for a review) have shown that internal two-body relaxation 
and tidal shocks cause stars to escape from clusters, even- 
tually leading to complete cluster dissolution. With the ex- 
ception of systems lying within a relatively narrow rang e 
of orbital parameters (see e.g. IVesperini fc Heggi3 1 19971 ). 
two-body relaxation is the dominant process driving the 
evaporation of stars from clusters. However the two-body 
relaxation timescale, Troiax, increases with both the mass 
(M) and the size (R) of the system, T^iax oc M^'^R^'^ 
(see e.g. iHeggie fc HutI |2003| ). For a cluster with an ini- 
tial mass of 10^ Mq, as adopted in the simulations pre- 
sented in Sect. |4l the mass loss induced by internal two- 
body relaxation in one Hubble time is negligible. For exam- 

! ile, the dissolution time pre dicted by A^-body simulations 
Baumgardt fc Makinol |2003| ) for a tidally limited cluster 
with initial mass moving on a circular orbit 

at 4 kpc from the Galactic center, is about 480 Gyr. Hence 
two-body relaxation cannot be the process responsible for 
the loss a significant fraction of the initial FG population, 
or for the evolution of the ratio of the number of SG to FG 
stars. 

A more promising mass reduction mechanism for the 
massive stellar systems considered in this study is the early 
mass loss associated with SN explosions in the FG popula- 
tion. In response to the loss of SN ejecta, a stellar system can 
expand beyond its tidal l imit (|Ch crnoff fc Wein berg 199^ ; 
iFukushige fc HeggiellT995h . losing stars from its outer lay- 
ers. The fraction of mass lost due to this expansion, and the 
implications for a cluster evolution, depend on the cluster 
structure, IMF and initial degree of mass segregation. 

For initially unsegregated, tidally limited clusters 
with a Salpeter IMF, it has been shown that the 
mass loss resulting from this mechanism is negligi- 
ble if the clusters are sufficiently centrally concen- 



trated, with concentrations (as measured by the con- 
centration paramete r c defined earlier) similar to thos e 
typical (1 < c < 2; [McLaughlin fc van der Marell (|2005|)) 
of Galactic globular cluste r s jV esperini fc Heggie 199^ 



McMiUan fc Portegies ZwartI [2003 : Baumg ardt fc Makinol 
2003D . However, low-concentration clusters can lose a sig- 



nificant fraction of their mass i n this way, a nd quickly dis- 
solve] Chcrnoff fc Shapiro 19871: IChernoff fc Weinberg 1990|; 
IFukushige fc Heggielll99a . see also IVesperini fc Zepj (I2003h 
for the possible implications of this dissolution on popula- 
tions of globular clusters.) 

Both the response of a cluster to this early mass loss 
and the critical concentration for its survival depend sensi- 
tively on the initial spatial distribution of massi ve stars. Sev- 
eral observational and theoretical studies (e.g. Hillenbrand! 
19971: iHillenbrand fc HartmannI Il998l; iFischer et al.l Il998l; 



de Griis et al' '2002"; 'Sirianni et al' '2002'; 'Goulierm is et al, 
2004; Stohc ct al. 2006; Bonnell ct al. 2001; McMilla n et al 



20071 ) have suggested that clusters may form with a signifi- 



cant degree of initial mass segregation imprinted by the star 
formation process, or creat ed dynamically very ear ly in their 
evolution. As discussed in IVesperini et al] l|2008al lbl). initial 
mass segregation has dramatic implications for early clus- 
ter evolution. Specifically, in an initially segregated cluster, 
the mass lost due to SN explosions is preferentially removed 
from the innermost regions of the cluster, and can lead to 
rapid and significant overall expansion, mass loss, and pos- 
sibly dissolution, even for an initially highly concentrated 
system. Our N-body simulations do not include primordial 
gas; the effect of primordial gas expulsion would add to that 
of SN ejecta in triggering the cluster early expa nsion (see 
e.g. iBoilv fc Kroupa 20031: iBaumgardt fc Kr oupa 2007) . 

In the following sections we explore the role of this early 
expansion on the evolution of multiple population clusters 
and, in particular, its implications for the evolution of the 
cluster stellar populations relative abundance. 

We begin with a general qualitative discussion of the 
scenario we have explored then illustrate the process using 
the results A^'-body simulations. 



6.1 Cluster mass loss 

In this section we outline the main elements of the dynamical 
mass loss scenario, the processes and timescales involved, 
and the effect on the numbers of stars in the two populations. 

The early evolution of the cluster is characterized by 
the rapid loss of FG stars. This is a consequence of clus- 
ter expansion in response to the dynamical heating from 
the loss of SN ejecta from the FG population: as the clus- 
ter expands beyond its tidal radius, the outer layers, which 
are predominantly populated by FG stars, are stripped. SG 
stars, which are mostly confined to the central regions, re- 
main approximately constant in number during this phase. 
As a result, the number ratio of SG to FG stars increases 
rapidly. Hereafter, we will refer to the number ratio of long- 
lived SG to FG main-sequence stars with masses in the range 
0.1 < M/Mq < 0.8 as /ms- 

The characteristic timescale for this early evolution, ex- 
pansion and stripping is the cluster dynamical time, Tdy„ oc 
although the precise duration of this phase may 
range from ~ 10 to lO'^ dynamical times depending on the 
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degree of initial mass segreRatio n and the amount of im pul- 
sive mass loss; see e.g. Fig. 3 in IVesperini et al] (|2008al 'l. 

The strength of this early expansion and the extent of 
the consequent loss of FG stars depend mainly on the degree 
of initial mass segregation, the amount of impulsive mass 
loss due SN and the structure of the cluster. As mentioned 
above, for clusters initially following a King density profile 
l|King II 19661 with concentrations similar to those currently 
observed in the Galactic globular cluster system initial mass 
segregation plays a key role in tri ggering a strong e xpansi on 
and early mass loss. As shown bv lVesperini et al] (|2008al lbh 
the heating due to a given amount of impulsive loss of SN 
ejecta is augmented by the preferential removal of the mass 
from the central regions. 

A second important initial parameter is the ratio of the 
cluster "King radius", ri<— that is the limiting radius of the 
King model used to describe the sytem- to the cluster's 
tidal radius, rtidai in the Galactic field: if rx/rtidai < 1, 
the cluster will not lose a significant number of stars until 
it fills its Roche lobe. For a cluster initially undefiUing their 
Roche lobe, /ms will start to evolve later and will eventually 
level off at a lower value than would be the case for a cluster 
with the same degree of initial mass segregation but initially 
filling the Roche lobe. 

As the early expansion phase slows and eventually stops 
so do the loss rate of FG stars and the time evolution of 
/ms. The subsequent "late" evolution is driven by two-body 
relaxation and proceeds on a timescale, Treiax , which is sig- 
nificantly longer than the dynamical timescale. Two-body 
relaxation is responsible for mixing the SG and FG popula- 
tions and for additional mass loss. 

For the massive clusters of principal interest in this pa- 
per, the total mass loss due to two-body relaxation is small; 
for example, u sing the mass loss rates deriv ed from N-body 
simulations bv lBaumgardt fc Makinol l|2003l ) and the current 
observational values of mass and galactocentric distance, we 
estimate that a cluster like lj Gen should have lost at most 
about 5-10 percent of its initial mass in one Hubble time by 
two-body relaxation. 

The study of less massive multiple population clusters 
undergoing a significant relaxation driven mass loss will be 
the subject of a future paper. We note, however, that star 
loss induced by two-body relaxation affects all cluster re- 
gions and not just the stars in the outer layers as in the early 
phases of the cluster evolution. Both SG and FG stars escape 
as a result of two-body relaxation with the result that the 
variation of the relative abundance of the two populations is 
significantly smaller during this late relaxation-driven phase 
than during the early dynamical cluster expansion. 

Although relaxation-driven stellar escape has a smaller 
effect on the relative abundance /ms of the two populations, 
the evolution of this ratio may be affected by mixing, an- 
other process driven by two-body relaxation. The relevant 
timescale in this case is the two-body relaxation time of 
the SG subsystem. As the low-mass SG stars evolve toward 
energy equipartition, the SG subsystem tends to expand, 
populating the outer regions of the cluster and mixing with 
FG stars. Once the two populations are mixed, even if the 
cluster is still undergoing a strong mass loss, the escaper 
population is no longer preferentially composed of FG stars 
and the number ratio of SG to FG stars stops increasing. 
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Figure 14. Initial density (normalized to the total central density 
po) profiles of FG (solid line) and SG (dashed line) stars. 

6.2 N-body simulations: initial conditions 

We now present the results of A'^-body simulations aimed at 
quantitatively exploring the scenario outlined above. 

We start our simulations with an initially segregated 
FG cluster of 25k stars. Mass segregation is set up by first 
letting the cluster evolve without including the effects of 
stellar evolution until the desired degree of segregation was 
reached by normal two-body relaxation. We emphasize that 
this is just the procedure used to generate a self-consistent 
initially segregated cluster model. It is not intended to model 
any stage of cluster evolution. The degree of initial segrega- 
tion is such that the ratio of the initial half-mass radius of 
the entire system {r^) to the half-mass radius of stars more 
massive than 1 Mq {rh,i) is rh/rhi — 1.5. 

Our model cluster is placed on a circular orbit at a 
galactocentric distance of Rg =4 kpc, with the Galac- 
tic tidal field modeled as a simple Keplerian potential of 
a point mass Mg ~ 4.3 x 10^° Mq (corresponding to a cir- 
cular velocity Vg ~ 220 km/s). We consider both a Roche- 
lobe-filing system {tk = ftidai), and three cases in which the 
system initially underfills its Roche lobe: rK/ftidai = 0.75, 
7'K/rtidai = 0.6 and rK/^tidai = 0.5. Hereafter we refer to 
these three simulations as SG-Rl, SG-R075, SG-R06, and 
SG-R05 respectively. 

SG stars are inserted into the simulation in a simplified 
way, by replacing, after about 30 Myr, the progenitor popu- 
lation of FG intermediate-mass stars (4 < M/Mq < 8) by a 
SG population having masses between 0.1 Mq and 0.8 Mq 
drawn from a Kroupa IMF. Following the results of the hy- 
drodynamical simulations presented in Sect.|4l newborn SG 
stars are concentrated in the innermost regions of the clus- 
ter. In the A''-body simulations, the SG stars are arbitrar- 
ily distributed according to a high-concentration Wq = 7 
(c = 1.5) King model, with a half-mass radius one-tenth 
of that of the initial FG system. The SG population is in- 
troduced with negligible velocity dispersion, and so under- 
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goes an initial rapid phase of violent relaxation. Largely for 
this reason, the details of initial SG density distribution are 
unimportant in determining the subsequent evolution of the 
SG/FG system. The adopted initial density profiles of the 
two stellar populations are shown in Fig. 1141 

As shown by an alytical calculat i ons an d A'^-body sim- 
ulations presented in lVesperini et al] l|2008al lbh. the heating 
due to mass loss from (FG) SN ejecta increases both with 
the degree of initial mass segregation of the objects losing 
the mass and with the amount of impulsive mass loss. In 
order to clarify the implications of a lower level of initial 
mass segregation of the FG population and/or of a smaller 
amount of early impulsive mass loss from FG SN ejecta, 
we have repeated the SG-Rl simulation after decreasing the 
velocities of all FG stars so to reduce the initial virial ratio 
by approximately 10 per cent and 30 per cent (changing di- 
rectly the initial virial ratio of the system allows us to easily 
focus on and explore the effect of a different initial degree of 
mass segregation while leaving all the other cluster proper- 
ties unchanged). Hereafter, we refer to these two simulations 
as SG-CIO and SG-C30. 

As described in the previous section, the scenario we 
are exploring is one in which cluster evolution is driven by 
multiple physical processes operating on different timescales. 
For the massive systems of interest here, all the relevant 
timescales (the dynamical time, the SG subsystem relax- 
ation time and the cluster relaxation time) are significantly 
different from one another. In particular, the early dynami- 
cal evolution phase of the cluster is quite distinct from the 
later phase of relaxation-driven mixing and mass loss. In or- 
der to maintain this clear separation between the early and 
late evolutionary phases for the relatively small-A system 
we can actually study, we have introduced into our A-body 
simulations a softening parameter equal to the average in- 
terparticle distance within the cluster half-mass radius in 
the calculation of the force. The effect of this softening is to 
increase the relaxation timescale, avoiding the overlap be- 
tween early rapid dynamical processes and late mixing and 
mass loss driven by two-body relaxation that would other- 
wise be found in the A-body models. 

All simulations w ere carried out u sing the 
Starlab package (|Portegies Zwart et al.l l200ll . 
http : //www ■ manybody . org/ ) accelerated by GRAPE-6 
special-purpose hardware (|Makino et al.l |2003| ). Starlab 
includes the effects of mass loss due to stellar evolution, 
two-body relaxation, and the influence of the Galactic tidal 
field. 

6.3 N-body simulations: results 

Fig. [T^] shows the time evolution of the number of SG, 
NsG,m3, and FG, Npcms main sequence stars with 0.1 < 
M/Mq < 0.8 for the SG-Rl simulation. The rapid prefer- 
ential loss of FG stars during the cluster early evolution, 
discussed in Sect. 16.11 is clearly evident. As the cluster ex- 
pands beyond its tidal radius in response to the loss of FG 
SN ejecta, its outer layers are stripped and a large fraction 
of the initial FG population is lost. The SG subsystem, on 
the other hand, is concentrated in the cluster inner regions 
and is largely unscathed by this early evolutionary process. 

The early cluster expansion and associated mass loss 
slow significantly as the cluster gradually enters the phase 
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Figure 15. Time evolution of the number of SG (dashed line), 
FG (solid line) main sequence stars with 0.1 < M/Mq < 0.8 for 
the SG-Rl simulation. 



driven by two-body relaxation. The time evolution of /ms = 
NsG,ms /NpG.ms for all the models explored in our simula- 
tions is shown in Fig. 1161 The early strong preferential loss of 
FG stars leads to a rapid increase of /ms; subsequently when 
the cluster expansion and mass loss slow down, the growth 
rate of /ms gradually decreases. Fig. [16] clearly illustrates 
the role of initial conditions in determining the current value 
of /ms and suggests a possible origin for cluster-to-cluster 
variation in the relative abundance of different populations. 
Our simulations indicate that the typical observed values of 
/ms ~ 0.5 — 1.5 are attained by clusters initially underfilling 
their Roche lobe and/or with moderate levels of initial mass 
segregation. 

Interestingly, our simulations also show possible evolu- 
tionary routes to the loss of most of the FG populations, 
leaving a cluster with a SG-dominated population. Initial 
conditions triggering the strong expansion necessary to strip 
most of the FG population are characterized by a larger de- 
gree of initial mass segregation and/or a larger amount of 
impulsive mass loss. The amount of impulsive mass loss will 
differ, for example, in clusters with the same stellar IMF 
and the same level of mass segregation but with different 
initial sizes (for example tidally truncated clusters at dif- 
ferent galactocentric distances) and therefore different dy- 
namical times. These clusters would lose the same amount 
of FG mass due to SN ejecta but for clusters with longer 
dynamical times, a larger fraction of this mass loss would 
occur in the impulsive regime, contributing to the heating 
that drives the cluster initial expansion. 

As the structure and the abundances of the FG and SG 
populations evolve, the relative spatial distribution of the 
two populations also undergo a significant variation. 

Fig. [T7] shows the evolution of the 10, 25, 50, 75 and 
90 per cent lagrangian radii for FG and SG main sequence 
stars with masses in the range 0.1 < M/Mq < 0.8 for the 
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Figure 16. Time evolution of /mSi the ratio of the number of 
SG to FG main sequence stars with 0.1 < M/Mq < 0.8, for all 
the Af-body simulations presented in the paper. 



SG-Rl simulation (repeated in this case without softening 
in order to properly interpret the evolutionary timescale in 
terms of the SG subsystem initial half-mass relaxation time). 

Fig. llSl shows the radial profile of /ms at three different 
times during the SG-Rl simulation. Our simulations predict 
that an initially strongly centrally peaked /ms (j) profile will 
evolve toward a profile that is fiat in the inner regions and 
decreasing in the outer parts. As the system evolves and the 
two populations continue to mix, the flat part of the pro- 
file (i.e. the region where the two populations have similar 
number density profiles) expands. 

For a detailed comparison between the level of mixing 
in real clusters and that predicted by our simulations, we 
need an estimate of the ratio of the cluster age, T, to the 
initial SG subsystem half-mass relaxation time, trh,SG(0). 
The lack of strong constraints on the initial conditions com- 
plicates such a comparison. However if we plausibly adopt 
the current cent ral relaxation time of uj Gen (~ 1 Gyr, see 
e.g. lHarris|[l99^ as an estimate of the initial half-mass re- 
laxation time of the SG subsystem, T /trh,SG(Q) ~ 10. The 
radial profile of /ms, predicted by our simulations for this 
value of T /trh,SG{0) is shown in Fig llSI It is characterized 
by a central flat region where SG and FG stars are mixed, 
and declines in the outer cluster regions still dominated by 
FG main sequence stars. 

We note that the presence of a declining outer part in 
the profile of fusir) is qualitatively consistent with the ob- 
served radial profile of t he number ratio of bl ue to red main 
sequence stars in uj Gen jSollima et al. II2007I ). Observations 
of the relative abundance of different populations in the in- 
ner regions of uj Gen and other multiple population clusters 
are needed to test the inner fiattening predicted by our sim- 
ulations. It is important to point out that while our pre- 
dictions concerning the shape and evolution of the /ms(7') 
profile are general, the quantitative characteristics of this 




d 



5 10 15 20 25 30 

t/trhSG(0) 

Figure 17. Time evolution of the 10, 25, 50, 75, 90 per cent 
lagrangian radii of the FG (solid line) and the SG (dashed line) 
population of main sequence stars with 0.1 < M/Mq < 0.8 for 
the SG-Rl simulation. trh,SG{0) is the initial half-mass relaxation 
time of the SG subsystem. 

profile (e.g. the slope of the outer parts and the extension 
of the fiat part of the profile) depend on several parameters 
including the initial relative concentration of the two popu- 
lations, the details of their initial density profiles, the initial 
degree of mass segregation, the number of stars left in the 
cluster after the early mass loss and the cluster relaxation 
time. We therefore emphasize that a quantitative compari- 
son between our theoretical /ms (r) profiles and observations 
must await the results of a more extended survey of simula- 
tions further exploring the space of initial parameters. 

Our simulations show that both /ms and the relative 
concentration of the two populations are significantly af- 
fected by dynamical evolution, and both depend on the clus- 
ter dynamical history and its initial conditions. We empha- 
size that the goal of the simulations reported here is to il- 
lustrate the key ingredients and dynamical processes behind 
the evolution of the relative number and spatial distribution 
of FG and SG stars. A more thorough survey of the param- 
eters relevant to the scenario presented here is needed for 
a detailed comparison with observational data, and is cur- 
rently in progress. The results will be presented elsewhere. 



7 CONCLUSIONS 

In this paper we have explored a scenario for the formation 
and the dynamical evolution of multiple stellar generations 
in globular clusters. A successful model for the origin and 
the evolution of clusters with multiple stellar generations 
needs to satisfy a number of observational constraints con- 
cerning the current chemical properties, relative abundance 
and structural properties of the different stellar populations. 

After reviewing the main observational results on mul- 
tiple population clusters, we have pointed out that the key 
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Figure 18. Radial profile of /ms at t/*rh,SG(0) = 0,10,25 
(dashed black line, dot-dashed blue line, solid red line) for the 
SG-Rl simulation. 

issues to address concern the general origin of SG stars 
(among which those with very high helium content observed 
in the most massive clusters are only a fraction) and the 
identification of the evolutionary processes leading to clus- 
ters with a similar number of SG and FG stars (or even to 
clusters dominated by SG stars). 

In order to address these different constraints and be 
able to model all the relevant aspects of multiple population 
clusters, we have combined the input from stellar evolution 
models to constrain the origin of the gas from which SG 
stars might form, hydrodynamical simulations to model the 
process of SG star formation and predict its initial structural 
properties and, finally, A'^-body simulations to explore the 
subsequent dynamical evolution of the two populations. 

The main steps and results of this investigation are as 
follows. 



1. Chemical Abundances 

We have reviewed all the issues concerning the super- AGB 
and massive AGB stars as the source of the polluted gas 
from which SG stars could form. The hypothesis that AGB 
stars might provide the gas for the SG formation has recently 
received some criticism in the literature based on possible 
inconsistencies between models and the observed chemical 
abundances in SG stars. However, in our review, we have 
pointed out that the results of chemical processing during 
the Hot Bottom Burning phases are highly afi'ected by the 
uncertainties in the basic input physics of these models, and 
that there is an ample range of plausible input parameters 
which can provide yields consistent with observations. In 
addition, recent models show that the super-AGB evolu- 
tion can produce the helium content required for the very 
helium-rich population mentioned above. During the super- 
AGB epoch, between the end of SN II explosions and the 



beginning of the SN la explosions, there are no mechanisms 
preventing the formation of a cooling flow and the forma- 
tions of the very helium-rich SG population from super-AGB 
winds. 



2. Hydrodynamical Simulations 

We have carried out a number of hydrodynamical simula- 
tions to explore the dynamics of gas from super-AGB and 
AGB stellar winds and the process of formation of the SG 
population in the assumption that the SNe II belonging to 
the FG population clear the GC of its pristine gas. 

(i) Our simulations show that this gas collects in a cooling 
flow into the innermost regions of the cluster where it forms 
SG stars. The cluster emerging from the hydrodynamical 
simulations is one with a SG strongly concentrated in the 
inner core of a more extended FG population. Initially the 
FG stars are the dominant stellar population with a total 
mass that must be about ten times larger than the total 
mass of the SG stars. The initial mass of the FG stars needs 
to be large enough to provide enough stellar mass return 
and form a substantial amount of SG stars. 

(ii) We have explored the dependence of our results on 
the initial mass of the cluster, the stellar IMF and the initial 
cluster size and concentration. SG stars are always concen- 
trated in the cluster innermost regions and the structure of 
the newly formed SG subsystem is largely independent of 
the cluster size and concentration. 

(iii) Starting with a initial FG mass ten times smaller 
than the value adopted in our standard model, leads to a 
SG total mass also ten times smaller but does not alter the 
shape of the SG density profile or its concentration. For sim- 
ulations with a flatter IMF the larger amount of gas supplied 
by AGB winds implies that a smaller star formation effi- 
ciency is needed to form the required mass of SG stars. If 
the requirement that the current mass of long-lived SG and 
FG stars are similar is maintained, a large FG initial mass 
is needed also for a fiatter IMF. In this case the number of 
long-lived FG stars is smaller, and a large value of Myg is 
still needed to have a current mass of long-lived FG stars 
consistent with observations. 

(iv) The effects of SNe la and additional energy sources 
on the SG formation process have been investigated. Our 
simulations show that, for a "normal" (e.g. Kroupa) IMF of 
the FG stars, the stellar explosions quickly evacuate most of 
the gas returned by the AGB stars, leading to a substantial 
drop in the gas content of the GC and to a suspension of 
the SG star formation. Given the uncertainties on the SN 
la rate, the interruption of the SF may occur in the time 
interval 40 < t < 100 Myr. Flat IMFs produce an increment 
in the amount of returned gas, leading to an enhancement of 
the radiative losses and, possibly, to the inability of the SNe 
la to get rid of the gas. In this case SG stars continue to form 
and their chemical abundances should show the signature of 
the pollution by gas from SN ejecta. This however does not 
seem to be confirmed observationally and flat IMFs are thus 
disfavoured. 

(v) We have further considered the possible presence of a 
diffuse energy source due to X-ray binaries and/or planetary 
nebulae. We show that for the standard model the cooling 
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flow and the associated SG star formation can be fialted by 
energy sources witli a total power larger than 10^^ erg s~^. 

(vi) We have also explored the possibility that the FG 
SNe II do not vent away all the pristine gas, and that part 
of it remains gravitationally bound to the cluster and falls 
back onto it after an initial ousting. The SG stars form- 
ing in the cluster central regions before this pristine gas 
reaches the cluster core originate from pure ejecta of the 
super- AGB stars which is very helium rich. These stars may 
form the very helium-rich population present in the most 
massive GCs. Once the pristine gas has fallen back in the 
cluster, SG stars form in an ambient medium in which the 
helium abundance is "diluted," and give rise to the moder- 
ately helium rich population. 

3. Stellar Dynamical Simulations 

The subsequent dynamical evolution of the cluster has been 
studied by means of A^'-body simulations. 

(i) One of the critical issues for viability of the scenario 
presented here is the identification of a mechanism leading 
to the escape of a large fraction of the initial FG population. 
For clusters as massive as required by our hydrodynamical 
models, the mass lost in one Hubble time due to two-body 
relaxation is negligible. 

Following observational and theoretical studies suggesting 
that initial mass segregation may be imprinted in clusters 
by the star formation process or produced dynamically early 
in the cluster evolution, we have studied the evolution of 
initially segregated clusters. In our simulations, a large frac- 
tion of FG cluster stars are lost early in the cluster evolution 
due to the expansion and stripping of the cluster outer lay- 
ers resulting from early mass loss associated with FG SN 
ejecta. Initial FG mass segregation plays a key role since 
the heating and the expansion due to the loss of SN ejecta 
are augmented by the preferential removal of the mass from 
the inner regions of the cluster. 

(ii) The population of escapers is dominated by FG stars. 
The SG population, initially concentrated in the innermost 
cluster regions, is largely unscathed by the early mass loss. 
We find that the early cluster evolution and mass loss can 
lead to a significant loss of FG stars and to current values 
of /ms, consistent with observations (0.5 < /ms < 1-5). We 
have also demonstrated possible evolutionary routes leading 
to the loss of most of the FG population and leaving a SG- 
dominated cluster. Clusters initially filling their Roche lobe 
with higher degrees of initial mass segregation and/or los- 
ing impulsively a larger fraction of SN ejecta undergo more 
expansion and may lose most of their FG stars. 

(iii) After the initial phase of expansion and mass loss, 
cluster evolution is driven by two-body relaxation which re- 
sults in the spatial mixing of the two populations. Following 
the hydrodynamical simulations, in the A-body initial con- 
ditions SG stars are concentrated in the inner regions of the 
cluster. The radial profile of the number ratio of SG to FG 
stars, fMs{r) is initially a decreasing function of radius. We 
find that, as the cluster evolves and the two populations 
mix, fMs{r) flattens in the inner parts of the cluster. The 
central flat region of the fMs{f') radial proflle expands on a 
relaxation timescale as mixing proceeds. In general, unless 
mixing is complete, we predict that /ms('') is characterized 



by a flat inner part and a declining portion in the outer 
regions of the cluster. 

A more extended survey of hydrodynamical and A-body 
simulations to further exlore the dynamics and the implica- 
tions of the scenario presented here as well as for a more de- 
tailed comparison with observations is currently in progress 
and will be presented in future papers. 
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APPENDIX 

The result of our simulations proved to be essentially inde- 
pendent of the value of v. To understand why, let us consider 
a one-zone toy model in which the gas density p and the stel- 
lar density p» evolve following the two equations 

-^-.p+s (13) 

dp* 2 

— ^up (14) 

The constant v regulates the SF efficiency and the constant 
s is a source term representing the rate at which the gas 
is replenished. The SF timescale t^i — l/fp depends on the 
gas density as the cooling time rather than the free fall time, 
which instead is the actual timescale encountered in our sim- 
ulations. With the present definition, however, it is possible 
to obtain an analytical solution to the above equations with- 
out invalidating our main conclusions (see Sect. [3l2|l . The 
solution of Eq. [TH] is 

P(i)= (^)° 'tanh(^) , (15) 
where r = (i/s)^"'^. Then the solution of Eg. [Til is 



..,)^(i)"'[(i)-..«..(i) 



(16) 



Defining the replenishing timescale as t^p = p/s, one obtains 
r = (tsitrp)'^'^ ■ In words, the growth timescale of the stars 
is given by the geometric mean of the SF and replenishing 
timescales, and is therefore of the same order of magnitude. 
This timescale is necessarily much shorter than the evolution 
time of the system over which egs. 1131 and 1141 are solved. In 
fact, these equations are meaningless for t/r < 1 because 
neither appreciable SF nor significant stellar mass loss can 
occur during a time interval that short. In the meaningful 
case in which t/r 3> 1, Eq. 1161 becomes 

p4t)^st. (17) 

The growth of the stellar population is thus independent of 
f and depends on the replenishing rate, a result confirmed 
by our simulations. 
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